******************* Table 2 Descriptive Statistics and Correlation******************
quietly zinb  multilateralagreement  new_DirectTransfer_dummy      YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy     YEAR1-YEAR6    )  forcevuong zip
asdoc sum multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity if e(sample), stat(N mean sd min max ) replace label
asdoc corr multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity if e(sample), label


******************* Table 3 Zero-Inflated Negative Binomial Models for ILAs******************
***** model 1-7 ****************
***** model 1 ZINB****************
zinb  multilateralagreement  new_DirectTransfer_dummy      YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy     YEAR1-YEAR6    )  forcevuong zip
est store a
***** model 2 ZINB****************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy     YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy  new_IndirectTransfer_dummy   YEAR1-YEAR6    )   forcevuong zip 
est store b
***** model 3 ZINB****************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust YEAR1-YEAR6    )   forcevuong zip
est store c
***** model 4 ZINB****************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
est store d
***** model 5 ZINB****************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    cluster(DYAD_id)   zip
est store e
esttab a b c d e using output_model1-5.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)  scalars(ll chi2  N_zero vuong  alpha)   replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  aic bic label
***** model 6 Hurdle NB****************
hnblogit  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 
est store f
esttab f using Table3_model6.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)  scalars(ll chi2 chi2_cp N_zero vuong)   replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  aic bic label
***** model 7 NB****************
nbreg  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , cluster(DYAD_id) 
est store G
esttab G using Table3_model7.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)  scalars(ll chi2 chi2_cp N_zero vuong)   replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  aic bic label
***** model 8 AME****************
zinb  multilateralagreement  i.new_DirectTransfer_dummy i.new_IndirectTransfer_dummy  i.had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( i.new_DirectTransfer_dummy i.new_IndirectTransfer_dummy   i.had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
eststo margin: margins, dydx(*) post
est store K
esttab K using Table3_model8.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)    replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  


///// Figure 3
gen IV1 = new_DirectTransfer_dummy
quietly zinb  multilateralagreement  IV1 new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( IV1 new_IndirectTransfer_dummy   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
margins, at(IV1=(0(1)1))  saving(file1, replace)
drop IV1
gen IV1 = new_IndirectTransfer_dummy
quietly zinb  multilateralagreement  new_DirectTransfer_dummy IV1  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy IV1   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
* same logit, different varnames
margins, at(IV1=(0(1)1)) saving(file2, replace)

grstyle init
grstyle set plain, grid
grstyle set color Set1
grstyle set legend 2
grstyle set graphsize 7cm 12cm
grstyle set size 10pt: heading
grstyle set size 8pt: subheading axis_title
grstyle set size 6pt: tick_label key_label
grstyle set symbolsize 3 4 5, pt
grstyle set linewidth .3pt: plineplot
grstyle set linewidth .4pt: legend axisline tick major_grid
grstyle set linewidth 0: pmark
grstyle set margin zero

combomarginsplot file1 file2, xtitle("")  ylabel( 0(1)3,format(%13.0fc) labsize(medlarge)) legend( size(medsmall) color(black) lcolor(white) fcolor(white)) graphregion(color(white)) ytitle("{stSerif:Predicted }" "{stSerif:Number of ILAs}", orientation(horizontal) size(large) margin(medium)) title("{bf:{stSerif:Figure 3.}} {stSerif:Effects of Direct and Indirect ANC on ILAs}", size(large) color(black) margin(medlarge))  plot1opt(symbol(S)  msize(medlarge) lpattern(dash)) legend(order(1 "Direct Agent Network" 2 "Indirect Agent Network") size(medsmall) color(black) lcolor(white) fcolor(white) margin(small)) xlabel( 0(1)1,format(%13.0fc) labsize(medlarge))

******************* The supplementary table A******************
********* Model A1. ZINB ***************************
zinb  multilateralagreement  i.new_DirectTransfer_dummy i.new_IndirectTransfer_dummy  ILA_lagged E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( i.new_DirectTransfer_dummy i.new_IndirectTransfer_dummy   ILA_lagged E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
est store S1
esttab  S1 using Supple_modelA1.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)    replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  
********* Model A2. ZINB ***************************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum wastewater_i industrial_wastegas_i  indussolidwaste_i industdust_i  pop_growth_i budget_growth_i gdp_growth_i cityi_size_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum wastewater_i industrial_wastegas_i indussolidwaste_i industdust_i pop_growth_i budget_growth_i gdp_growth_i cityi_size_km2  proximity   YEAR1-YEAR6)    forcevuong zip
est store S2
********* Model A3. AME ***************************
eststo margin: margins, dydx(*) post
est store S3
********* Model A4. ZINB ***************************
zinb  multilateralagreement  new_DirectTransfer_dummy new_IndirectTransfer_dummy  had_been_cooper_dum wastewater_j industrial_wastegas_j  indussolidwaste_j industdust_j  pop_growth_j budget_growth_j gdp_growth_j cityj_size_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy new_IndirectTransfer_dummy   had_been_cooper_dum wastewater_j industrial_wastegas_j indussolidwaste_j industdust_j pop_growth_j budget_growth_j gdp_growth_j cityj_size_km2  proximity   YEAR1-YEAR6)    forcevuong zip
est store S4
********* Model A5. AME ***************************
eststo margin: margins, dydx(*) post
est store S5
esttab S2 S3 S4 S5 using Supple_modelA2-5.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)    replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  
********* Model A6. ZINB ***************************
zinb  multilateralagreement  i.new_DirectTransfer_dummy##i.new_IndirectTransfer_dummy had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( i.new_DirectTransfer_dummy##i.new_IndirectTransfer_dummy  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6    )    forcevuong zip
est store S6
esttab  S6 using Supple_modelA6.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)    replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  

********* Model A7. ZINB ***************************
sort DYAD_id year
bysort DYAD_id: gen new_DirectTransfer_dummy_lag = new_DirectTransfer_dummy[_n-1]
bysort DYAD_id: gen new_IndirectTransfer_dummy_lag = new_IndirectTransfer_dummy[_n-1]
zinb  multilateralagreement  new_DirectTransfer_dummy_lag new_IndirectTransfer_dummy_lag  had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust  abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity  YEAR1-YEAR6 , inflate( new_DirectTransfer_dummy_lag new_IndirectTransfer_dummy_lag   had_been_cooper_dum E_abs_WasteWater E_abs_WasteGas_万亿立方英尺 E_abs_SolidWast_百万吨 E_abs_indusdust abs_pop_growth_hundred abs_budget_growth_hundred abs_gdp_growth_hundred abs_citysize_km2  proximity   YEAR1-YEAR6    )    forcevuong zip
est store S7
esttab  S7 using Supple_modelA7.rtf ,  rtf star(* 0.1 ** 0.05 *** 0.01)    replace  nogaps drop (YEAR1 YEAR2 YEAR3 YEAR4 YEAR5 YEAR6)  se(%8.3f)  b(%8.3f)  

